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The second (unphysical) critical charge in the 3-body quantum Coulomb system 
of a nucleus of positive charge Z and mass nip, and two electrons, predicted by 
F Stillinger has been calculated to be equal to = 0.904854 and Z^^ = 0.905138 
for infinite and hnite (proton) mass mp, respectively. It is shown that in both cases, 
the ground state energy E(Z) (analytically continued beyond the first critical charge 
Zc, for which the ionization energy vanishes, to ReZ < Zc) has a square-root branch 
point with exponent 3/2 at Z = Zb in the complex Z-plane. Based on analytic 
continuation, the second, excited, spin-singlet bound state of negative hydrogen ion 
H“ is predicted to be at -0.51554 a.u. (-0.51531 a.u. for the finite proton mass nip). 
The first critical charge Zc is found accurately for a finite proton mass mp in the 
Lagrange mesh method, = 0.911069 724655. 
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We consider the 3-body quantum Coulomb system of two electrons (e, m) and a (heavy) 
positive charge Z of mass rup, {Z,mp), where m <C mp. It is one of the most fundamental 
systems in theoretical physics. In atomic units |e| = 1. Usually, this system is called the 
helium isoelectronic sequence, among many different names used. We denote this system 
as {Z, e, e) and prefer to call it the two-electron sequence. This system is fundamental: it 
includes the species like H~, He, Li+ etc, which play exceptionally important role in Nature. 
The non-relativistic Hamiltonian has the form 
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when written in atomic units with electronic mass m = 1, where ri(r 2 ) is the distance 
from the charge Z to the hrst (second) electron, ri 2 is interelectron distance, Ai(A 2 ) is 
the Laplacian which describes the kinetic energy of the hrst (second) electron, and A is 
the Laplacian which describes the kinetic energy of charge Z. The conhguration space of 
o is R®, but after separation of centre-of-mass motion, we arrive at a six-dimensional 
conhguration space of the relative motion. Note that in the static limit mp —)■ oo the 
non-relativistic Hamiltonian ([T|) degenerates to the form 


H = -5(A. 
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Due to charge quantization, the physical system occurs when Z takes integer values, Z = 
1,2,.... At Z > 1 the energy of a bound state E{Z) is a smooth function proportional 


Q 


to (—Z^) at large Z. Long ago it was given a rigorous mathematical proof in [lH that for 
rrtp = 00 the function E(Z')IZ‘^ is the analytic function in vicinity of 1/Z = 0. There is 
still unresolved challenge to hnd the radius of convergence in 1/Z expansion and establish 
nature of the singularity(ies) on its circle of convergence. The lower edge of continuum 
(the threshold energy) is given by the ground state energy of the Z-hydrogen atom (Z, e), 
Eth = ~^^ 2(l+mp) . It is well-known that there exists a certain charge Z - called the critical 
charge Zc - for which the ground state energy of a two-electron sequence degenerates with 
the lower edge of continuum, thus, spontaneous ionization becomes possible. Quite recently, 
in the static limit ([2]), the critical charge was accurately calculated (with triple basis sets 


containing up to 2276 terms, see [2j) and verihed in ind^endent calculation (using the 


Lagrange mesh method with lattice size 92 x 92 x 20, see S) within 12 signihcant hgures 


zoo 


0.911028 224 077 . 


( 3 ) 
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Q 


It was proved in that at Z = the Schrodinger equation has normalizable eigenfunction. 
In this result was verified - it was explicitly stated that the trial function (triple basis 
set containing 2276 terms) is normalizable and localized, < ri >~ l/ai ~ 1 a.u. and 
< r 2 >~ l/a 2 ~ 5a.u. , see below (6)-(7). It implies that the probability to meet electron 
at inhnity is equal to zero, thus, the spontaneous ionization is absent. 


The critical charge for the case of the finite mass rup was studied in for some 
particular values of m/mp and in 6| for > 1 and m = 1. However, only recently, it was 
calculated accurately for the mass of proton, nip = 1836.15267245me, where rUe = 1 was 
set for the electron mass, using a 2856-term trial function (the convergence in energy was 
checked up to 11 s.d.) 


0 


= 0.911069 7(3) , (4) 

which now we conhrm (and improve) in the Lagrange mesh method (see for details 3) 
(the convergence in energy is checked up to 12 s.d. including, the maximal lattice size is 
92 X 92 X 20) 

Z^^ = 0.911069 724 655, (5) 

Note that finite-mass effects are sufficiently small ~ 10“^ changing the hfth hgure in the 
critical charge, see (|3]) and ([5]). 

It is evident that when Z gets slightly smaller than Z^. the bound (stationary) state 
becomes quasi-stationary state (or Gamow state), characterized by complex energy, for 
discussion see [w, 11|. Its wavefunction is complex: it is mostly localized in the Coulomb well 
and it corresponds to outgoing spherical wave at large r 2 (ri) at hxed ri(r 2 ), ~ 

The imaginary part of energy is defined by tunneling rate from the Coulomb well to infinity. 
In principle, it can be found via multidimensional WKB method, it should be exponentially 

b 

small, ~ (Zc — Z)“e (Zc-z)=^ where 6, c > 0 at Z —)■ Z^—. We are not aware about any 
concrete theoretical calculations done so far. Hence, we have two different spectral problems 
which can be summarized as follows: (i) for Z > Zc we have a problem of bound states 
looking in the Schrodinger equation for solutions in the Hilbert space with (real) energies 
E = Ey{Z), (ii) for Z < Zc we have a problem of quasi-stationary states looking in the 
Schrodinger equation for solutions corresponding to outgoing spherical wave with (complex) 
energies E = E^{Z). Thus, the boundary conditions for these two problems are completely 
different. In principle, there is no reason that these two analytic functions should coincide 
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E^{Z) = E^{Z). Furthermore, the function E^{Z) approaching to Z —>■ Z^— is complex 
with exponentially small imaginary part (it indicates to the essential singularity at Z = Z^, 
while the Ey{Z) is real and approaching from the right to Z —)■ Zc-\- behaves linearly j4 . 
It looks like a type of Kosterlitz-Thouless phase transition (of the inhnite order) occurs at 
Z = Zc- the functions E^{Z) and Ey{Z) and all their derivatives match at Z = Z^- The 
goal of our study is to hnd the function Ey{Z) - the energy of the ground state - its analytic 
properties. We will not touch the function E^{Z) either theoretically, or numeric ally . 

In remarkable papers of 1966 and 1974, F.H. Stillinger and D.K. Stillinger 12|, ll3|, made 
an extended analysis of the ground state energy obtained variationally as a function of 
continuous charge Z and predicted the existence of the second critical charge Zb - the 
charge for which the ground state eigenfunction looses its normalizability, see 12| . The trial 
function that was used for the analysis, was the celebrated Hylleraas-Eckart-Chandrasekhar 
function with a property of clusterization - in average one electron is closer to the charge 
center than the other, they are called innermost and outermost electrons. 


'^triai{ri,r2,ri2) = To(ri,r2) + To(r2,ri) , 4'o(n,?^2) = e ai ^ a2 , (6) 

where ai^ 2 {Z) are two variational parameters, both real, found in the minimization of the 
energy functional for a hxed Z. The Z dependence of the parameters ai ^2 was quite non¬ 
trivial. It must be noted the important technical moment which made the analysis feasible: 
all integrals involved are evaluated analytically and are rational functions in ai^ 2 - Hence, 
the variational energy is a simple analytic function of parameters ai ^2 being the ratio of two 
polynomials in variables ai^ 2 - It allows analytic considerations of the variational energy^. 
Stillinger and Stillinger showed the existence of two distinct critical charges: (i) Zc{= 0.9538) 
for which the ionization energy vanishes, hence, the ground state energy is equal to —Z^/2 
and coincides (degenerates) with the lower bound of continuum (the ground state energy 
of the Zc-Hydrogen atom, (Zc,e)), and (ii) Zb{= 0.9276) for which the trial function ([6]) 
becomes non-normalizable: minimal ai (or 0 : 2 ) vanishes, for illustration see Fig. [T] It was 
also demonstrated that the variational ground state energy E^ar = E{Z) is a regular function 
at Z = Zc, while at Z = it has a square-root branch point with exponent 3/2 with the 
branch cut going along real axis to the negative direction, Z G (— 00 , Zb]- Recently their 


^ The same is correct for the function 0, see below 
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analysis was extended: it was checked that, assuming the parameters ai ^2 in (E]) can be 
complex, the minimum of energy functional is always reached for real ai 2 for Z > Zq. 

The ground state variational energy can be constructed for a more general trial function 
with the electronic correlation explicitly included |l^ . 


d^o(n,r2,ri2) = (l + cri2)e , ai ^ a2 , 


(7) 


(c.f. (E])), where ai^ 2 {Z),'^{Z),c{Z) are four variational parameters. Similarly to ([ 6 ]), the 
ground state was characterized by two distinct critical charges, Zc 7 ^ Zb, where Zc{= 0.9195) 
is regular point while Zb{=^ .8684) is the square-root branch point with exponent 3/2. Thus, 
the introduction of two extra parameters c, 7 to the trial function ([ 6 ]) shifts significantly the 
position of each critical point without changing the analytic properties. Note that both 
trial functions (Ej), ([7]) lead to sufficiently accurate ground energy for Z = 1 (which is a 
weakly-bound state): —0.5133 and —0.5260, respectively, while E^xact = —0.52775a.u. and 
for Z = 2 : —2.8757 and —2.9019, respectively, while Eexact = —2.9037 a.u. In fact, (E]) and 
(17|) are the most accurate trial functions among two and four-parametric trial functions for 
any integer Z G [1, 00 ): they give 1-2 and 2-3 correct decimal digits in ground state energy, 
respectively. Thus, both variational energies provide highly accurate, uniform(!) approxi¬ 
mation of the exact energy in Z G [1, cxo). Striking common property of both approximate 
variational energies as analytic functions of Z is the absence of singularity at Z = Zc. Note 
that the existence of two critical charges was also demonstrated in the Born-Oppenheimer 
approximation rrip —?■ 0 (two-center, H^-type case) jl5 |. 

The analysis of both variational energies Eyar{Z) obtained with (E]), ([7]) did not indicate 
the existence of other square-root branch points at finite Z on the complex Z-plane jl^ . 
Hence, one can draw the conclusion, in accordance with the Landau-Zener theory of the 
level crossings (for discussion see e.g. Q, Q), that the point Zb seems to be a natural 
candidate for the point of the level crossing of the two lowest spin-singlet, l^S* and 2^5* states. 
However, the analytic continuation of the energy to the real positive axis of the second sheet 
of the Riemann surface in Z in the vicinity of Zb is not fully supported by the analytic 
continuation of trial functions (E]), ([7]). Although these functions remain normalizable after 
analytic continuation, they do not develop a nodal surface, thus, they are NON-orthogonal 


to the ones corresponding to t 
Extending the analysis of 


le grou nd state at given Z and need to be orthogonalized. 


12 


, [l^, one can derive that for both cases ([ 6 ]) and ([7]) the 
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variational energy in Z has two essentially different functional expansions near and Zg, 
namely: 

E{Z) = - Y + - ^c) + a2(Z - Z,f + a3(Z - + ... , (8) 

which is the Taylor expansion: it indicates the absence of a singularity at Z = Zc, and 

E{Z) = bo + bi{Z — Zb) + ci{Z — Zb)^ + b2{Z — ZbY + C2{Z — Zb)^ + ■ ■ ■ , (9) 


for Z > Zb, which is the so-called Puiseux expansion (the expansion in fractional degrees). 
The latter indicates the existence of the square-root branch point with exponent 3/2. In 
absence of other singularities in the vicinity of Z^ other than at Z = Zb, the radius of 
convergence of the expansion ([H]) has to be equal to (Zc — Zb)- Since the continuum starts 
at Eth = —Zl/2, the presence of the bound state at E indicates the phenomenon 


of the existence of a bound state embedded to continuum 


18| 


It seems true inside the 


accuracy of a given variational procedure. It must be noted that even highly accurate, 
normalizable trial function found at Z > Zc does not guarantee the existence of the solution 
of the Schrodinger equation in the Hilbert space for Z < Z^ after analytic continuation. 
Thus, this state is not necessarily bound. 

It is worth emphasizing that the expansions (IH]) and ([9]) are exact functionally since they 
are derived from variational energies emerging from ([ 6 ]) and ([^ after the exact minimization: 
minimizing-the-energy variational parameters are found in a form of expansions near Zc and 
near Zb+, respectively. 

It seems natural to check validity of the expansions (|H]) and (E]) for more accurate energies 
than ones found using (E]), ((Tj). In order to do it we calculate E{Z) in 12 points at Z G [0.91,1] 
and Z = 2 in the Lagrange mesh method accurate to not less than 9 s.d. (see Table [I] for 
illustration). Then taking into account the points from the domain Z G 1] only, 

we make a fit using a general terminated Puiseux expansion, 

N 


E — cin{Z — Zb 


kn/2+e„ 


n=0 


for different N assuming that exponents are growing with u, searching for optimal values of 
the parameters an, e„. As for initial values of parameters, those we use, they are ones taken 


^ Recently, this unexpected statement was supported in highly accurate variational study [2| : for Z slightly 
smaller that Zc, Z < Zc normalizable trial functions guarantee convergence in energy up to 14 decimal 
figures without loosing normalizability! 
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z 

E (a.u.) 

fit (ground state) 

fit (2nd branch) 

2.00 

-2.903 724 4 

-2.903 719 

-2.201 927 

1.00 

-0.527 751 0 

-0.527 749 

-0.515 540 

0.94 

-0.449 669 0 

-0.449 668 

-0.447 012 

0.92 

-0.425 485 3 

-0.425 485 

-0.424 741 

ZEBMD (ggg p) 

-0.414 986 2 

-0.414 986 

-0.414 793 

0.91 

-0.413 799 2 

-0.413 799 

-0.413 652 


Table I: Ground state energy E for a two-electron system at rup = oo for selected values of Z 
found in the Lagrange mesh method and rounded to 7 d.d., compared with the fit (jlOp : the result 
of the analytic continuation of (I1U|) to the second branch is shown in the 3rd column. 


from the expansion ([9]) for ([6|), Finally, we arrive at 


E[(f{Z) = -0.407924 - 1.12347 (Z-Z^) - 0.197785 (Z-Zb)^/^ _ 0.752842 (Z-Z^)^ 


0.108259 (Z - Zb)®/2 _ 0.014135 (Z-Zb)^+ 0.00854 (Z 
+ 0.00483 (Z - Zs)^ - 0.000056 (Z - Zb)9/2 , 


ZbY^^ 


( 10 ) 


with the critical charge 

Zb = 0.904854 , (11) 

where as the resnlt of the fit e„ tnrned to be of order 10 “® — 10 “^ as well as the parameter 
oi. It tnrns ont that is artifact of htting procednre: all these parameters e^, oi can be placed 
eqnal to zero withont redncing the qnality of the fit! The expression flTOl) reprodnces np to 
a portion of 10“® — 10“^ in energies at 12 points in Z more or less eqnally distribnted in 
Z G [Zc, 1.] and also the point Z = 0.91 < Zc (!) , for illnstration see Table I (and Table II, 
see below). The qnality of the fit drops with the increase of Z, see Table I. The nnmerical 
coefficients in the expansion decrease with the increase of the order of terms. It may be 
considered as the indication to the finite radins of convergence of corresponding nnmerical 
series of coefficients and thns the Pnisenx expansion ([9]) itself. We mnst emphasize that 
the critical charge Z^ is nnphysical: it corresponds to analytically-continned gronnd state 
energy to the domain of qnasi-stationary states! 

In the recent analysis of 1/Z expansion 19j it was shown that the first ten decimal 


digits in all the first 401 coefficients are fonnd correctly in 


20 j: their weighted snms 












Ep{Z) = reproduce subsequently the first ten decimal digits in ground state 

energy of two-electron sequence at Z = 1,2,... 10, found in 2]J]. Checking convergence 
of 1/Z expansion we calculate weighted sums Ep for Z < 1, see Table Ull and compare 
them with the result of fit (ITOll and with energies obtained in Lagrange mesh calculations 
on the lattice 92 x 92 x 20 . Comparison of Ep and indicates to a striking agreement 

between them: they differ in ~ 2 ■ 10 “^ for Z = 0 . 95 , in ~ 6 ■ 10 “® for Z = Z^ ^ 0.911 and 
~ 4 • 10 “^ for Z = 0 . 905 . Hence, the fit ffTOj) is in agreement with perturbation theory in 
1/Z. In a similar way Ep agrees with accurate calculations of the ground state energy E in 
the Lagrange mesh method even for Z < Zc- 


z 

Ep 

E[{f (unD 

E (Lagr.mesh) 

0.905 

-0.408 045 

-0.408 089 

-0.408 li 

0.908 

-0.411 478 

-0.411 501 

-0.411 502f 

0.909 

-0.412 629 

-0.412 648 

-0.412 648 5f 

0.91 

-0.413 783 

-0.413 799 

-0.413 799 2 

z- 

-0.414 972 

-0.414 986 

-0.414 986 2 

0.92 

-0.425 482 

-0.425 485 

-0.425 485 3 

0.93 

-0.437 450 

-0.437 451 

-0.437 451 3 

0.94 

-0.449 669 

-0.449 669 

-0.449 669 0 

0.95 

-0.462 125 

-0.462 125 

-0.462 124 7 


Table II: nip = cx3 case: Weighted sum Ep{Z) (a.u.) with 1/Z coefficients 6^ (first column), see 
text; Energies (1101) in a.u. (second column) from the fit based on a terminated Puiseux 

expansion (IinD; Energies E (a.u.) in Lagrange mesh method (rounded to 7 s.d., third column), the 
results marked by f contain confident decimal figures only 


Moreover, it has to be mentioned a remarkable fact that expanding flT0|) a.t Z = oo the 
established coefficients in 1/Z- expansion Cio-^o ai'e reproduced with relative accuracy 10“^ 
(equivalently, the first two significant digits) (9|. 

Similar analysis to (fTOj) can be made for the hnite-mass case nip. Making a fit for the 
energies at Z E [Z/^^, 1] (see Table HTT]) we arrive at 

EUit .rrxp) ^ -0.408019 -1.123511 (Z - Zp) - 0.198005 (Z - Zpf^'^ - 0.751842 (Z - Zpf 
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z 

E (a.u.) 

fit (ground state) 

fit (2nd branch) 

2.00 

-2.903 304 6 

-2.903 178 

-2.196 303 

1.00 

-0.527 445 9 

-0.527 444 

-0.515 312 

0.96 

-0.474 536 5 

-0.474 537 

-0.469 305 

0.92 

-0.425 242 4 

-0.425 242 

-0.424 519 

0.911069725 (see ©) 

-0.414 798 1 

-0.414 798 

-0.414 617 

0.91 

-0.413 564 0 

-0.413 564 

-0.413 430 


Table III: Ground state energy E for a two-electron system at nip for selected values of Z, found 
in the Lagrange mesh method, rounded to 7 d.d. and compared with the fit (I12p : the result of 
analytic continuation of (I12|) to the second branch is shown in 3rd column. 

- 0.101848 {Z - - 0.02171 {Z - Zef + 0.0039 {Z - Zsfl'^ + 0.0126 {Z - Zef 

- 0.0028 (Z - ZbT'^ , (12) 

with the critical charge 

Zb = 0.90514 , (13) 

(cf ([HD). Note that the parameters of the fit ffTU]) obtained in the static limit change very 
little in going from infinite to finite proton mass, see flT^ . Both critical charges ■^c, Zb 
increase slightly for the case of the finite proton mass. 

It can be immediately seen from Tables I. lIIIl that based on the analytic continuation in Z 
of the ht (ITOH . flT^ (it corresponds to the change of sign in front of the terms of half-integer 
degrees), the existence of the second bound state 2^5' of negative hydrogen ion is predicted, 

= -0.51554 a.M. , = -0.51531 a.u. , (14) 

for the both infinite and finite proton mass, respectively. Transition energies are 

= -0.0122 a.M. , ^^ 5^215 = -0.0121 a.M., (15) 

respectively. H“ is a ’’maximally” strongly-correlated Coulomb system, thus, it is an example 
of a general 3-body Coulomb problem, it should be treated in full generality. So far, we are 
unable to find accurately the wavefunction for the 2^5* state. It is evident it should not be of 
the type (ls2s) as for He, Li’*' and other two-electron positive ions. Likely, relativistic effects 
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Figure 1: Z-complex plane of the ground state energy E = Ey, see text, the branch cut shown 
by bold wide black line; vertical line passing through Zc indicates the line of discontinuity (non- 
analyticity) between E^{Z) and E^{Z) (see text). 


(spin-orbit, spin-spin interactions, radiative effects) to energy of the 
should not be more than ~ 10“‘^a.u. similarly to ones for Helium 22 1 


owest states 


Note that conceptually the prediction flT^ is in agreement with D.R. Yafaev’s rigorous 


mathematical statement 


23| about the hnite number of bound states at Z = 1 but in 


contradiction to the widely-known theoretical statement by R.N. Hill 2J, [2^ about the 
existence of a single, spin-singlet bound state of H“ only. Mathematical justihcation of 
the latter result is absent (to the best of the present author’s knowledge). It seems it is 
a challenge to check the prediction (IT^ experimentally since the transition —?■ 2^5* is 

the forbidden type transition, hence, it should be a type of multi-photon transition and/or 
bound-free-bound transitions, which is difficult to observe. In principle, the spectra of H“ 
contains one more, but metastable spin-triplet 2^P state 2^. Note that for helium atom, 
the main modes are (permitted) electric-dipole transitions 2^5* —?■ 2^P and 2^P —)■ HS", see 
271. |28| . Experimental study of the (forbidden) transition HR —)■ 2^5* for helium atom was 
done only recently using the Doppler-free two-photon spectroscopy, see [29|. 

In present paper we studied the domain of bound states of three-body Coulomb problem 
Z > Zf. and the analytic continuation to the complex Z-plane of the ground state l^R, 
in particular, to Z < Zc. We showed that the second critical charge predicted in [l2| is 
unphysical: it appears in analytic continuation of the ground state energy to Z < Zc, thus, 
it can not be measured. It seems important to carry out a similar study for the domain of 
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quasi-stationary states Z < Zc, in particular, checking the analytic discontinuity Z = 
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